Breakup of 8 B and the Sn astrophysical factor revisited 
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Existing experimental data for the breakup of 8 B at energies from 30 to 1000 MeV/nucleon on 
light through heavy targets are analyzed in detail in terms of an extended Glauber model. The 
predictions of the model are in excellent agreement with independent reaction data (reaction cross 
sections and parallel momentum distributions for core like fragments). Final state interactions 
have been included in the Coulomb dissociation component. We extract asymptotic normalization 
coefficients (ANC) from which the astrophysical factor £17(0) for the key reaction for solar neutrino 
production, 7 Be(p,7) 8 B, can be evaluated. Glauber model calculations using different effective 
interactions give consistent, though slightly different results. The differences give a measure of the 
precision one can expect from the method. The unweighted average of all ANCs extracted leads to 
517(0) = 18.7 ± 1.9 eVb. The results of this new analysis are compared with the earlier one. They 
are consistent with the values from most direct measurements and other indirect methods. 
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I. INTRODUCTION 

The major source of the high-energy neutrinos ob- 
served by the solar neutrino detectors is 8 B, produced 
in the 7 Be(p,7) 8 B reaction [l( at the end of the pp III 
chain. The recent results from Superkamiokande 2] and 
SNO shift the interest for a precise determination of 
the rate of this reaction from the problem of the exis- 
tence of the solar neutrino deficit and of the neutrino 
oscillations to that of putting stringent constraints on 
the different scenarios that explain them. There were 
many recent determinations of Sit, but its precise value 
is still controversial. In particular, there is a discrepancy 
between the value found in one direct measurement and 
most of the results from indirect measurements. 

Recently we have proposed an indirect method to ex- 
tract astrophysical S factors from one-nucleon-removal 
(or breakup) reactions of loosely bound nuclei at inter- 
mediate energies 0, 0- It is based on the recognition 
that the structure of halo nuclei is dominated by one 
or two nucleons orbiting a core (see for example @, 
and references therein). Consequently, we use the fact 
that the breakup of halo or loosely bound nuclei is es- 
sentially a peripheral process, and therefore the breakup 
cross-sections can give information about the wave func- 
tion of the last proton at large distances from the core. 
More precisely, we determine asymptotic normalization 
coefficients (ANCs) from a comparison of the experimen- 
tal data with calculations. Then, these ANCs are suffi- 
cient to determine the astrophysical S factors for radia- 
tive proton capture reactions. The approach offers an 
alternative and complementary technique to extracting 
ANCs from transfer reactions 8] , an alternative particu- 
larly well adapted to rare isotope beams produced using 
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fragmentation. 

In this paper we discuss the use of existing experi- 
mental data on 8 B breaku p a t energies between 30 and 
1000 MeV/nucleon 0, [l(J lH [H [13 to determine the 
astrophysical factor Sn- The calculations presented in 
on this subject were extended and refined. First, the 
Coulomb part of the dissociation cross section was mod- 
ified by including the final state interaction into the cal- 
culations. Second, new data on the breakup of 8 B are 
analyzed [ill [I2> E3] ■ Third and most important, a new 
set of calculations for the breakup of 8 B were made us- 
ing different effective nucleon-nucleon (NN) interactions. 
Each of the new effective interactions considered, which 
do not involve any new parameters, give consistent re- 
sults for all experiments, but the average ANCs found 
are slightly different from one interaction to another. We 
interpret these differences as a measure of the accuracy of 
the present (and possibly other) indirect method(s). Fi- 
nally, a brief comparison with results of direct measure- 
ments and of other determinations of Sn using indirect 
methods is made. 



II. FROM BREAKUP CROSS SECTIONS TO 
ANCS 

In the breakup (one-nucleon removal reactions) of 
loosely bound nuclei at intermediate energies, a nucleus 
B = (Ap), where B is a bound state of the core A and 
the nucleon p, is produced by fragmentation from a pri- 
mary beam, separated and then used to bombard a sec- 
ondary target. In inclusive measurements, the core A is 
detected, measuring its parallel and transverse momenta 
and eventually the gamma-rays emitted from its deexci- 
tation. Spectroscopic information can be extracted from 
these experiments, such as the orbital momentum of the 
relative motion of the nucleon and the contribution of 
different orbitals (from the momentum distributions) and 
the contribution of different core states (from the coin- 
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cidences with gamma-rays). Typically the experimen- 
tal results are compared with calculations using Glauber 
models. The integrated cross sections have been used to 
extract absolute spectroscopic factors Q or the ANC . 
We have shown that the latter approach has the advan- 
tage that it is independent of the geometry of the proton 
binding potential, an important feature for exotic nuclei 
for which the geometry of the mean field is not neces- 
sarily well known. The ANC for the nuclear system 
A + p «-> B specifies the amplitude of the tail of the 
overlap function of the bound state B in the two-body 
channel (Ap) §]. This ANC is enough to determine the 
direct (non-resonant) contribution to the astrophysical S 
factor for the radiative proton capture reaction A(p, j)B 
which is a highly peripheral process due to the Coulomb 
barrier and the low energies in the entrance channel. Us- 
ing this strategy we described the breakup of 8 B in terms 
of an extended Glauber model. The 8 B projectile (made 
of a proton and the 7 Be core) is moving on a straight 
line trajectory and each part is interacting independently 
with the target. The breakup cross sections depend on 
the proton-target and core-target interactions and on the 
relative p-core motion. 

The wave function of the ground state of 8 B is a mix- 
ture of lj»3/2 and lpi/2 orbitals, around a 7 Be core: 

\ s B(g. S .) > =A P3/2 [ 7 Be(3/2-)®p 3/2 )] 2+ + 

A Pl/2 [ 7 Be(3/2~)® Pl/2 )] 2+ + 

A e [ 7 Be*(l/2~)®p 3/2 )} 2+ +... (1) 

where Ai are the spectroscopic amplitudes of the various 
components. The first two terms represent the proton 
in the I-P3/2 and lpx/2 orbitals, respectively, coupled to 
the ground state of 7 Be. The third term corresponds 
to the proton being coupled to the first excited state of 
the 7 Be core, at E* =0.429 MeV. Basic shell model argu- 
ments suggest that the IP3/2 term dominates, and only 
a small lpi/2 admixture exists. Recently, in the study 
of its mirror nucleus 8 Li, we disentangled for the first 
time these two contributions and found their ratio to be 

Kj A l 3/ , = C lJ C l^ = °- 13 (2) M Only these 
two terms contribute in the radiative capture process. 
However, all three terms contribute in the breakup pro- 
cess, with the third one identified in 8 B breakup through 
coincidences with gamma-rays |13|. It does not con- 
tribute in the radiative capture, but its contribution has 
to be evaluated and subtracted from all the other in- 
clusive breakup data. From the breakup cross section 
to the excited state in 7 Be, a(exc.) = 12(3) mb and 
a (tot) — 94(9) mb measured at 936 MeV/nucleon, we 
found Cg/C t 2 ot = 0.16(4), a value consistent with that 
found in the original analysis in |13| and which, subse- 
quently, was used to correct for the contribution of core 
excitation in all other breakup data analyzed here. These 
two findings together establish the wave function of the 
ground state of 8 B, up to an overall multiplicative factor. 

The calculated one-proton removal cross sections and 
the momentum distributions are given by the incoherent 



superposition of the single particle contributions from the 
different parts of the wave function weighted by the re- 
spective spectroscopic factors 

= X! S ( c ' nl 3) a s P (nlj)- (2) 

In inclusive measurements, such as those analyzed here, 
the proton is not detected, therefore the calculated cross 
sections u sp (nlj) contain a stripping term (the loosely 
bound proton is absorbed by the target and the 7 Be core 
is scattered and detected) , a diffraction dissociation term 
(the proton is scattered away by the target, the 7 Be core 
is scattered by the target and is detected) and a Coulomb 
dissociation term [l6j 

/>oo 

<r sp = / 2nbdb(P st r(b) + P d iff(b)) + a Cou i (3) 
Jo 

These terms were calculated using the extended 
Glauber model detailed elsewhere 0, 03 • S-matrix el- 
ements have been calculated in the eikonal approxima- 
tion up to the second order (T^j to assure convergence. 
This convergence was checked with calculations for other 
quantities, for example proton-target reaction cross sec- 
tions as a function of energy, and compared with data 
available from literature I20l. 



III. RESULTS USING DIFFERENT NN 
INTERACTIONS 

In calculations we assume a structure of the projec- 
tile given by Eq. ^ with the spectroscopic factors, or 
the ANCs, to be determined from the comparison of the 
measured cross sections (from which the contribution of 
the 7 Be core excitation was removed as described above) 
with those calculated as an incoherent superposition of 
single particle cross sections 

cr_ip = (S P3/2 + S Pl/2 )a sp = (Cl 3/2 + Cl i/2 )(j sp /b 2 p , (4) 

where b n ij are the asymptotic normalization coefficients 
of the normalized single particle radial wave functions 
f n ij(r) calculated in a spherical Woods-Saxon potential 
of a given geometry and with the depth adjusted to re- 
produce the experimental proton binding energy of 8 B, 
S p = 0.137 MeV. They are essentially equal for the IP3/2 
and lpi/2 orbitals (b p ), as are the single particle breakup 
cross sections o~ sp . The sum of the spectroscopic factors 
or the sum of the asymptotic normalization coefficients 
Cf ot = Cp + C| can thus be extracted by compar- 
ing the experimental one-proton removal cross sections 
with the calculations. The 8 B ANC, C t 2 ot , is extracted 
from existing breakup data at energies between 30-1000 
MeV/nucleon and on different targets ranging from C to 
Pb [I ESIH IH El- Figure ffla) shows the one-proton 
removal cross sections for various targets and incident 
energies. One can notice the large range of cross sections 
and the variation with the energy for different targets. 
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FIG. 1: a) The cross sections determined from the 
breakup of 8 B at 30-1000 MeV/u on C, Al, Sn and Pb tar- 
gets at various energies [9-13] used in this study. b)The 
ANCs determined from the breakup of 8 B using the JLM 
effective interaction. The error bars of the individual 
points contain the experimental and theoretical uncer- 
tainties. The dashed line shows the average and the 
hatched area is the standard deviation. 



Two approaches were used to evaluate the S'-matrices 
needed in the calculations. The first is a potential ap- 
proach. To obtain the folded potentials for the proton- 
target and core-target interactions we used the effective 
nucleon-nucleon interaction of Jcukenne, Lejeune and 
Mahaux (JLM) [2lJ and Hartree-Fock-Bogolyubov den- 
sities carefully adjusted to correctly reproduce the exper- 
imental binding energy of each nucleus. In an extensive 
study of the elastic scattering of loosely bound p-shell 
nuclei around 10 MeV/u [22 , we found that renormal- 
ized double folded potentials with this effective interac- 



tion provide a good description of the data. We found 
there that a large renormalization is needed for the real 
part of the potential, but no renormalization is needed for 
the imaginary part of the potential. In the present calcu- 
lations we assume that no renormalization of the imag- 
inary part is needed at all energies. We used the JLM 
interaction for energies below 285 MeV/nucleon only. 

Before comparing the experimental and calculated in- 
tegrated cross sections, we checked that we can repro- 
duce all other available experimental observables with 
our model. This was crucial before proceeding with the 
calculations. In Figure [21 we show that parallel momen- 
tum distributions measured at 41 MeV/nucleon on one 
low Z (Be) and one high Z (Au) target |23j and on the 12 C 
target at 936 MeV/nucleon (calculated with appropriate 
technique for high energy, as discussed below) for both 
the ground state and excited state components are well 
reproduced. Similarly, the transverse momentum distri- 
butions are well reproduced. The model also reproduces 
well the relative fraction of stripping/diffraction dissoci- 
ation disentangled first by Negoita et al. on Si targets 
at 28-38 MeV/nucleon (as can be seen in Fig. 6 of that 
reference) and more recently by Enders et al. on C at 
76 MeV/nucleon For the latter the calculations give 
&str = 80 mb, Udiff+c = 50 mb, to compare with the ex- 
perimental results <T s tr — 93(16) mb, Odiff+c — 37(13) 
mb. In Figure |3| we show, for the case of the breakup 
of 8 B on C targets that the reaction is essentially pe- 
ripheral. The stripping and the nuclear diffraction disso- 
ciation probabilities as a function of the proton impact 
parameter, s, are calculated at four energies. While these 
probabilities are peaked outside the radius of the 7 Be core 
(vertical line) in all cases, it is clear that the interior con- 
tributes and should be carefully considered. The figure 
also shows the variation with energy of the relative im- 
portance of the two nuclear mechanisms: the diffraction 
dissociation (dashed line) is dominant at lower energies 
and its role decreases with increasing energy where strip- 
ping (full line) becomes dominant. The comparison of 
the results of the present calculations with the results 
of the simpler black disk model shows that the interior 
plays the crucial role in describing correctly the wings 
of the parallel momentum distributions (see Fig. 3 of 
Ref. PI). An analysis like the one presented in Fig. 2 of 
Ref. [J shows that there is an energy window E/A=25- 
150 MeV/nucleon for which the breakup of 8 B is mostly 
peripheral even on the lightest targets. For the heavier 
targets this is always the case, due to the dominance of 
the Coulomb component. 

The data considered were taken on C targets at 76 
[HI, 142, 285 [13 and 936 A MeV (exps. no. 1-4, 
in order), on Al at 285 A MeV (U3, exp. #5), on Si at 
28^ 35, 38 (H, exps. #6,8,9), 29*, 39, 49 and 56 A MeV 
(Ha, exp. #7,10-12), on Sn at 142 and 285 A MeV (H3, 
exps. #13,14) and on Pb targets at 142 and 285 A MeV 
(E3, exps. #15,16). 

From the analysis with the JLM interaction of all 
experiments up to 285 A MeV we find ANCs consis- 
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have, therefore, performed calculations using bare ampli- 
tudes resulting from perturbation theory [25|], as well as 
renormalized E2 and El amplitudes. No significant dif- 
ferences were found in the extracted ANCs with these 
two versions, and the values reported here are those ob- 
tained without any renormalization. The Coulomb term 
in the breakup cross section is particularly important for 
heavy targets where it becomes dominant. The value 
found above for the ANC is in very good agreement with 
that determined before using the peripheral proton trans- 
fer reactions 10 B( 7 Be, 8 B) 9 Be and 14 N( 7 Be, 8 B) 13 C at 12 
MeV/nucleon ^ Cf ot (p) = 0.449 ±0.045 fni" 1 and with 
that obtained from the study of the mirror neutron trans- 
fer reaction ( 7 Li, 8 Li) C? ot {n) = 0.455 ± 0.047 fin" 1 0. 
They agree very well, in spite of the differences in the 
energy ranges and in the reaction mechanisms involved. 
The ANC extracted with JLM leads to the astrophysical 
factor Sn(Q) = 17.5 ± 1.8 eV- b for the key reaction for 
solar neutrino production 7 Be(p,7) 8 B. 
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FIG. 2: The parallel momentum distributions deter- 
mined from the breakup of 8 B on Be and Au targets at 41 
MeV/u :23] and on C at 936 MeV/u for the both g.s. 
and core excitation components. Final state interaction 
is included in the Coulomb calculations. The total (full 
lines) and the components are shown: stripping (dashed) 
and diffraction (dotted), Coulomb (dash-dotted), or as 
labelled on each curve. 



tent with a constant value (Fig. ^ b) with an average 
C? ot (JLM) = 0.454 ± 0.048 fin -1 . Compared with Ref. 
P| , we include the newer measurements by 0, 0] . An- 
other distinction is that we have included the final state 
interaction in the calculation of the Coulomb dissociation 
component of the one-proton removal cross section. El 
and E2 amplitudes have been included as in the earlier 
calculation, except that distorted waves, not plane waves, 
were taken in the p+ 7 Be final channel for the calculation 
of the matrix elements. The distorted waves were cal- 
culated numerically in the same potential that was used 
to bind the proton p around the 7 Be core in the ground 
state of 8 B. Differences occur between the calculated am- 
plitudes with the two approaches especially for low rela- 
tive momenta, but their influence on the final integrated 
result is relatively small due to the extra q 2 factor that 
weights their contribution to the integrated cross section. 
However, the inclusion of distorted waves increases the 
asymmetry in the parallel momentum distribution due 
to an increased E1-E2 interference effect as can be seen 
in the upper right panel in Fig|3 It has been suggested 
|24| that asymmetries observed in the fragment parallel 
momentum distributions in the Coulomb dissociation of 
8 B on heavy targets could be reproduced with an overall 
renormalization of 1.22 and of 0.7 for the E2 matrix el- 
ements calculated in first order perturbation theory. We 
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FIG. 3: The breakup probability profiles as a function 
of the impact paramenter s for the breakup of 8 B on 
C targets at four different energies. The stripping (full 
lines) and the diffraction dissociation (dashed lines) com- 
ponents are shown. The vertical line shows the position 
of the 7 Be core rms radius. 

In a second approach, the Glauber model in the opti- 
cal limit |23 was used. The breakup process is treated 
as multiple elementary interactions between the partners' 
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nucleons. The total NN cross sections and the scattering 
amplitudes are taken from literature. Calculations were 
done for all the experiments in the energy range 50-1000 
MeV/nucleon using a constant ("standard") finite range 
of 1.5 fm, as well as sp ecific ranges in each NN channel as 
suggested by Ray [2I|. No new parameters were adjusted. 
For details on the procedure see [f|- For all the effective 



1j 




10 12 14 16 

experiment 

FIG. 4: The ANCs determined from the breakup of 8 B 
at 30-1000 MeV/u on C, Al, Sn and Pb targets, using 
three NN effective interactions: JLM (squares), "stan- 
dard" (circles) and " Ray" (triangles). See text for details. 
The dashed, dotted and dash-dotted lines are the average 
of the JLM, "standard" and "Ray" values, respectively. 



taken in various laboratories, at different energies, with 
varying methods and the calculations also used different 
techniques. Therefore, we believe that the results form a 
statistical ensemble with many and randomly occurring 
error sources, for which the average and the standard de- 
viation around the average give a reasonable description 
of the ANC and its error. 

In Ref. [3l| the authors study the same data of 8 B 
breakup on the C target and find a larger value for the 
ANC than the one wc published previously in p|. They 
use a different strategy for the calculations where they as- 
sume a wave function for the 8 B g.s. from nuclear struc- 
ture calculations and a geometry of the proton binding 
potential that they do not question. Then, the compar- 
ison with the experiment gives them a quenching factor 
R s of unexplained origin in that paper (but of great sig- 
nificance if its connection with short range correlations 
inside nuclei is confirmed) . On the other hand they com- 
pare their result for one single target with the full aver- 
age from our calculations. A direct comparison with the 
individual ANCs or with the average of our results for 
the breakup on the C target only (available in Table I 
of our Ref. Q) would have led to agreement. Later [ll| 
they find full agreement with us Q for the breakup of 
9 C where we use essentially the same techniques. Also, 
our examination of different theoretical reaction mod- 
els above indicates that a quenching factor R s = 0.88 
may not be precise enough to consider it different from 
unity. A recent study of 23 cases of one-neutron removal 
cross sections at similar energies 17] found no quenching 
R ave = 0.98 ±0.16. 



NN interactions we checked that they correctly describe 
complementary data, like proton-target and 7 Be-target 
elastic and total reaction cross sections, where available. 
Understandably, calculations with zero-range and with 
2.5 fm range give too small or too large cross sections, re- 
spectively, and were not retained. Data at energies higher 
than 50 A MeV were selected. We did not include the 
measurements of Ref. [28j at 1440 MeV/nucleon and of 
Ref. [2|j at 1471 MeV/nucleon (highest energy points in 
Fig. njt) , because at those very large energies the breakup 
is no longer peripheral and the extraction of an ANC 
may not be the most appropriate. However, the results 
obtained from the analysis of these two cases are fully 
consistent with the others. 

For each of the two NN interactions wc find that all ex- 
periments give consistent ANCs (Fig. , but the average 
values obtained are slightly different: Cf ot (" standard" ) = 
0.503±0.032fm" 1 and C? ot ("Ray") = 0.517±0.041 fm" 1 . 
These differ by 11% and 13% respectively from the JLM 
value. We find no argument to determine which value is 
best. If we take the unweighted average of all 31 determi- 
nations we find an ANC C t 2 ot (awe) = 0.483 ± 0.050 far 1 
that leads to 5i 7 (0) = 18.7±1.9 eV- b. The uncertainties 
quoted are only the standard deviation of the individual 
values around the averages, with no experimental errors 
included. The experimental data considered here were 



IV. CONCLUSIONS 



In conclusion, we show that the breakup of 8 B at in- 
termediate energies can be used to obtain the Sn astro- 
physical factor at stellar energies. Very difficult direct 
measurements are complemented by reactions using sec- 
ondary beams of exotic nuclei obtained from fragmenta- 
tion and seeking the relevant ANCs, rather than a com- 
plete knowledge of the ground state wave function of 8 B. 
In addition, the indirect ANC method is subject to dif- 
ferent systematic errors than direct measurements. 

There were many recent determinations of this key as- 
trophysical factor Sn, but its precise value is still con- 
troversial. Our result is in agreement with those from all 
indirect methods and with most of the direct determi- 
nations (see the discussions in [32, IH, 03 El] ) > but one 
which stands out in its claim of a larger value and very 
small error .36]. The value obtained as an average of all 
ANCs found in the present study S 17 (0) = 18.7 ± 1.9 
eV b, is virtually equal with the most probable values 
extracted in Ref. [H S 17 (0) = 18.6 ± 1.2 (stat) ±1.0 
(theor) eV b and in Ref. S 17 (0) = 18.6 ± 0.4 (stat) 
±1.1 (theor) eV b from statistical analyses of all mu- 
tually consistent results, including the reanalysis of data 
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from direct measurements |3jj with a different extrapola- 
tion at low energies. Our results from the use of different 
NN interactions reminds us of the fact that the precision 
of indirect methods depends not only on the precision of 
the experiments but also on the accuracy of the calcula- 
tions. These findings may give a measure of the present 
status for break-up reactions, indicating that accuracies 
to +/-10% can be obtained. 
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